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1  ACCOMPLISHMENTS:  EXECUTIVE  SUMMARY  OF  RESEARCH 

1.1  Technical  Objective  and  Expected  Payoff 

The  enemy  threat  to  GPS  is  evolving  in  two  areas:  enemy  use  of  GPS  and  the  enemy’s  ability  to  deny 
GPS  to  the  United  States.  Worldwide  military  use  of  GPS  is  evolving  due  to  the  wide  availability 
of  commercial  GPS  receivers,  and  the  widespread  knowledge  of  the  force  enhancement  capabilities 
offered  by  GPS.  The  jamming  threat  is  a  concern  because  of  the  physical  design  of  the  GPS  system. 

The  received  power  from  the  GPS  satellites  is  approximately  -157  dBW.  Many  jammers  available  on 
the  arms  market  today  either  already  eover  the  GPS  frequencies,  or  can  be  modified  to  do  so. 

The  technical  objective  is  to  identify,  evaluate,  develop,  and  demonstrate  innovative  anti-jam  (AJ) 
processors  that  are  effective  against  intentional  and  unintentional  interference  to  GPS  user  equipment. 

This  entailed  a  system-level  study  of  the  trade-offs  between  protection  performance  versus  hardware 
and  computational  complexity  for  various  integrated  levels  of  processing.  The  primary  end  objective 
is  the  development  of  fully  digital  AJ  solutions  that  are  integrated  into  the  GPS  receiver  to  achieve 
higher,  yet  affordable  AJ  capabilities  for  GPS  user  equipment  on  navy  cruisers,  aircraft,  and  airborne 
weapons  systems .  However,  AJ  preprocessors  which  may  be  implemented  as  retro-fits  to  either  existing 
or  currently  planned  GPS  receivers  have  been  developed  as  well. 

The  expected  payoff  arises  from  the  proven  ability  of  GPS  to  provide  significant  force  enhancement 
capability.  This  force  enhancement  capability  has  been  demonstrated  in  every  U.S.  military  operation 
since  and  including  the  Gulf  War.  The  success  of  GPS  in  these  operations  has  led.  to  an  increasing 
Department  of  Defense  (DoD)  and  Allied  reliance  on  GPS.  To  date,  over  200,000  GPS  receivers  have 
been  delivered  to  the  DoD  and  its  Allies.  GPS  uses  include  navigation  targeting,  sensor  aiding,  weapons 
delivery,  timing,  and  battlefield  command/control.  A  key  performance  feature  of  the  algorithms  that 
have  been  developed  is  the  ability  to  defeat  electronic  attack,  including  the  smart  or  spoofing  jammer. 

The  time  is  right  for  a  fully  digital  solution  so  that  the  expected  pay-off  should  be 
very  substantial.  Prom  a  signalling  viewpoint,  GPS  is  a  Code  Division  Multiple  Access  (CDMA) 
system  wherein  each  satellite  is  assigned  a  different  code.  The  Third  Generation  (3G)  of  commercial 
cellular  CDMA  systems  is  currently  under  development.  To  support  wireless  access  to  the  Internet,, 
for  the  downloading  of  video  in  real-time,  for  example,  it  is  planned  that  3G  CDMA  systems  will 
accommodate  bandwidths  up  to  20  MHz  -  same  as  the  bandwidth  occupied  by  the  P  code  signal! 

Texas  Instruments  (TI)  and  other  companies  are  now  producing  DSP  chips  that  do  CDMA  processing 
commensurate  with  this  bandwidth.  In  addition,  commercial  CDMA  systems  involve  correlations  with 
a  so-called  “long”  code  that  is  similar  in  nature  (roughly. speaking)  to  the  P  code.  Further,  part  of 
3G  CDMA  will  operate  in  a  band  centered  at  1.9  GHz,  not  very  far  in  frequency  from  LI. 

The  point  is  that  in  the  very  near  future  fully  digital  solutions  effecting  AJ  for  GPS  may  be  feasibly 
implemented  on  commercially  available  DSP  chips  developed  for  3G  CDMA.1  These  DSP  chips  will  be 
mass  produced  facilitating  a  cost  effective  solution. 

The  overall  goal  of  this  project  is  to  develop  an  affordable  GPS  system-level  solution  that  protects 
DoD’s  and  U.S.  Allies’  continuous  operational  use  of  GPS. 

2  Military  Problem/Deficiency  Addressed  and  Resulting  Capability 
Enhancements 

2.1  Deficiency  Addressed 

The  success  of  GPS  in  Desert  Storm  and  other  operations  has  increased  DoD’s  commitment  to  equip¬ 
ping  weapon  systems  with  GPS  receivers.  However,  this  commitment  has  been  accompanied  by  a 

*It  is  interesting  to  note  that  just  as  GPS  chip  sets  are  coming  down  in  price  enough  to  be  feasible  for  use  in  cell 
phones  for  911  Emergency  Geolocation,  the  DSP  chips  embedded  in  cell  phones  are  becoming  advanced  enough  to  process 
signals  at  the  P  code  rate. 
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heightened  awareness  of  the  potential  vulnerabilities  of  GPS  to  both  intentional  and  unintentional 
RF  interference.  Limited  Electronic  Counter  Counter  Measures  have  been  developed  over  the  past 
two  decades  primarily  based  on  the  null  steering  Controlled  Reception  Pattern  Antenna  (CRPA)  and 
associated  Antenna  Electronics  (AE).  However,  this  spatial  .filter  technology  is  limited  in  its  ability 
to  cancel  jammers  due  to  component,  device,  and  technology  limitations,  as  well  as  the  increasing 
potential  jamming  threat. 

For  example,  the  CRPA  AE  system  was  one  of  the  earliest  anti-jam  spatial  filters  for  GPS,  a  seven 
element  null  steering  antenna  array.  The  principle  used  in  the  CRPA  AE  system  (deployed  on  F-16C 
fighters)  to  cancel  interference  is  power  minimization.  Power  minimization  is  premised  on  the  fact 
that  prior  to  correlation,  the  C/A  code  signal  from  any  one  GPS  satellite  is  on  the  order  of  30  dB 
below  the  noise  floor.  In  power  minimization  based  null  steering,  the  weight  for  the  reference  antenna 
is  set  to  unity  and  the  remaining  weights  applied  to  the  auxiliary  antennas  are  found  as  those  which 
drive  the  output  power  of  the  beamformeras  close  to  the  noise  floor  as  possible.  This  adaptively  steers 
nulls  in. the  directions  of  those  interferes  whose  power  levels  are  above  the  noise  floor. 

The  CRPA  AE  system  employed  vector  modulators  to  effect  the  weights  applied  to  the  auxiliary 
antennas  at  RF.  Vector  modulators  are  digitally  controlled  analog  devices  that  adjust  the  in-phase 
and  quadrature  output  of  an  antenna  at  RF.  One  of  the  primary  problems  with  the  CRPA  AE  system 
is  that  it  used  analog  power  measurements  to  drive  the  weight  adaptations;  this  has  problems  dealing 
with  pulsed  jammers,  among  other  things.  The  CRPA  AE  system  also  used  a  crude  coordinate 
descent  algorithm  in  adapting  the  weights,  which  converged  very  slowly  and  sometimes  exhibited, 
erratic  behavior. 

The  GAS-1  AE  system  represented  a  substantial  improvement  over  the  CRPA  AE  system.  The 
GAS-1  AE  system  was  partly  digital;  the  weight  updates  were  computed  digitally  using  correlations 
measured  from  sampled  data.  However,  similar  to  the  CRPA  AE  system,  the  GAS-1  AE  system 
employed  vector  modulators  to  adjust  the  weights  of  the  auxiliary  antennas  at  RF.  Vector  modulators 
do  not  allow  precise  control  of  the  weight  values,  thereby  impeding  convergence  and  limiting  null 
depths.  In  addition,  the  system  is  a  null  steering  antenna  array,  and  thereby  has  limited  degrees  of 
freedom  for  canceling  interferers,  especially  since  the  number  of  antennas  is  only  seven. 

2.1.1  Shortcomings  of  Null  Steering  Antennas 

A  limitation  of  null  steering  antenna  technology  is  that  space-only  processing  requires  the  placing 
of  a  spatial  null  in  the  direction  of  each  narrowband  interferer  as  well  as  each  wideband  interferer. 
This  leads  to  two  problems.  First,  the  maximum  number  of  interferers  that  can  be  spatially  nulled 
is  M  —  1,  where  M  is  the  number  of  antennas.  Of  course,  given  that  the  array  aperture  in  the  GPS 
application  is  typically  on  the  order  of  a  wavelength,  two  interferers  may  be  relatively  widely-spaced 
in  terms  of  physical  angle  but  closely-spaced  in  terms  of  beamwidths.  A  single  null  may  therefore 
potentially  take  out  two  interferers  (or  more  depending  on  their  locations).  Notwithstanding,  in  the 
GPS  application,  the  number  of  antennas  is  small  due  to  cost  considerations,  size  limitations,  power 
consumption,  etc.  Thus,  cancellation  of  interferers  solely  through  spatial  nulling  consumes  precious 
degrees  of  freedom  and  thereby  has  limited  AJ  capability. 

The  second  problem  is  that  if  an  interferer  and  a  particular  GPS  satellite  are  closely-spaced  in 
angle,  the  formation  of  a  spatial  null  towards  that  interferer  may  drop  the  gain  in  the  direction  of 
that  GPS  satellite  so  low  that  it  is  rendered  useless.  Again,  due  to  the  aforementioned  small  aperture 
problem,  an  interferer  and  GPS  satellite  direction  may  be  relatively  widely-spaced  in  terms  of  physical 
angle  but  closely-spaced  in  terms  of  beamwidths.  In  order  for  the  GPS  receiver  to  provide  accurate 
navigation  information,  it  is  necessary  to  track  the  signals  from  at  least  four  different  GPS  satellites. 
In  order  to  reduce  errors  from  various  sources,  it  is  desirable  to  track  the  respective  signals  from  a 
larger  number  of  GPS  satellites.  Thus,  if  certain  types  of  interferers  can  be  canceled  by  means  other 
than  spatial  null  steering,  this  would  be  highly  advantageous. 
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2.1.2  Inherent  Limitation  of  Power  Minimization  Based  AJ 

Again,  both  the  CRPA  AE  and  the  GAS-1  AE  system  employ  the  principle  of  power  minimization 
which  attempts  to  drive  the  output  power  of  the  beamformer  as  close  to  the  noise  floor  as  possible.  In 
so  doing,  nulls  are  only  steered  in  the  directions  of  those  interferers  whose  power  levels  are  above  the 
noise  floor.  This  is  another  key  shortcoming  of  current  technology.  Smart  jammers  (delay  spoofers, 
for  example)  can  lie  below  the  noise  floor  and  yet  totally  disrupt  subsequent  tracking  and  timing 
calculations  related  to  the  GPS  signals.  The  AJ  schemes  presented  for  this  effort  incorporate  additional 
layers  of  advanced  signal  processing,  over  and  above  power  minimization,  based  on  the  known  Gold 
Code,  for  example,  in  order  to  cancel  jammers  below  the  noise  floor. 

This  work  is  critical  since  it  is  expected  that  the  number  of  potential  military  and  commercial 
interferers  will  increase  along  with  more  and  new  jamming  threats  that  will  likely  be  fielded  in  response 
to  the  U.S.  fielding  of  numerous  GPS  equipped  weapon  systems.  Thus,  innovative  solutions  are 
needed  for  canceling  larger  numbers  of  jammers/interferences  while  minimizing  the  impact  on  the  GPS 
signal  as  best  as  possible.  Fully  digital  solutions,  ideally  integrated  into  the  GPS  receiver  itself,  offer 
greater  flexibility,  greater  accuracy,  and  a  larger  number  of  degrees  of  freedom  for  AJ  and  interference 
suppression. 

2.2  Resulting  Capability  Enhancements 

The  AJ  processors  developed  through  this  effort  are  capable  of  suppressing-  a  mixture  of  interferer 
types  including  multiple  broadband  Gaussian  noise  interferers.  This  type  of  interferer.  is  a  significant 
threat  as  it  is  possible  to  distribute  the  jammer  power  across  the  20.46  MHz  bandwidth  of  the  standard 
military  P(Y)  code  signal  centered  at  the  GPS  LI  frequency,  for  example,  at  a  low  enough  level  that  it 
is  difficult  to  detect  and  identify,  and  yet  the  noise  background  is  increased  enough  to  disrupt  tracking 
of  the  GPS  signal.  This  is  particularly  true  if  the  enemy  distributes  (spatially)  a  number  of  relatively 
low  power  broadband  jammers  over  the  combat  arena. 

The  processor  also  provides  suppression  capabilities  against  other  types  of  interferers,  both  inten¬ 
tional  and  inadvertent.  These  include  continuous  wave  (CW),  swept  CW,  pulsed  CW,  phase  shift 
keying  (PSK),  pseudo-noise  signals  (20  MHZ  bandwidth),  and  narrowband  and  wideband  frequency 
modulated  signals.  The  interferers  may  be  located  anywhere  within  or  adjacent  to  the  20.46  MHz 
bandwidths  centered  at  the  GPS  LI  frequency  of  1575.42  MHz  or  the  L2  frequency  of  1227.60  MHz, 
and  distributed  anywhere  over  2n  steradians  of  solid  angle  centered  at  zenith  relative  to  the  local 
horizontal  plane  of  a  GPS  receiving  antenna. 

In  addition  to  providing  interference  suppression,  the  processor  allows  reception  of  GPS  satellite 
signals  in  a  stressed  environment  by  maximizing  the  signal  power  to  interference-plus-noise  power 
ratio  (SINR)  for  acquisition  and  tracking  of  the  GPS  signal.  Various  user  environments  have  been 
addressed.  A  top-down  requirements  study  of  the  number  of  antennas  and  number  of  taps  per  antenna 
has  been  conducted  under  various  user  environments  ranging  from  a  foot  soldier  in  the  field,  to  a  navy 
cruiser,  to  a  high  performance  fighter  aircraft.  Various  levels  of  integrated  processing  have  been 
developed  thereby  offering  an  array  of  solutions  that  may  be  assessed  through  cost  versus  performance 
considerations. 

3  Integrated  AJ  Processor  Realizations 

We  have  developed  four  different  realizations  of  an  AJ  processor  for  GPS.  The  listing  below  is  in  order 
of  increasing  complexity  for  increased  protection  capability. 

•  Stand-Alone  Power  Minimization  Space-Time  Preprocessor  (Section  3.3) 

•  Per  Satellite  Power  Minimization  Space-Time  Preprocessor  (Section  3.5) 
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•  Multiple  Power  Minimization  Space-Time  Preprocessors  Followed  by  Per  Satellite  SINR  Maxi¬ 
mization  Combining  (Section  3.6) 

•  Power  Minimization  Space-Time  Preprocessor  Followed  by  Per  Satellite  MMSE  Receiver  Based 
on  Known  C/A  Code  (Section  3.7) 

The  first  AJ  processor  allows  current  GPS  receivers  to  be  simply  retro-fitted  with  a  front-end 
device.  As  a  practical  matter,  though,  in  order  to  truly  avoid  any  modification  to  the  GPS  receiver 
itself,  the  output  of  the  space-time  pre-processor  must  undergo  D/A  conversion  followed  by  mixing 
back  up  to  RF.  This  is  obviously  counter-productive  when  employing  a  digital  solution  but  necessary 
in  cases  where  a  simple  retro-fit  is  the  mandated  solution. 

The  other  three  A J  processors  generate  a  different  space-time  filtered  output  for  each  GPS  satellite 
in  the  field  of  view  (FOV).  Each  of  these  three  AJ  processors  may  be  implemented  in  a  stand-alone 
fashion  by  summing  the  outputs  for  the  different  satellites  ahd  feeding  the  resulting  sum  signal  to 
the  GPS  receiver.  However,  they  are  more  optimally  implemented  as  an  integrated  component  of  the 
overall  GPS  receiver  wherein  the  space-time  filtered  output  for  the  fc-th  GPS  satellite  is  fed  directly 
to  the  correlators  and  attendant  timing  operations  matched  to  the  fc-th  satellite. 

Each  of  the  three  more  advanced  AJ  processors  requires  an  estimate  of  the  Direction-of- Arrival 
(DOA)  vector  for  the  k- th  GPS  satellite  in  the  FOV.  The  m-th  element  of  the  M  x  1  DOA  vector 
for  a  given  GPS  satellite  is  the  gain  and  phase  of  the  GPS  signal  received  from  that  satellite  at  the 
m-th  antenna,  relative  to  the  gain  and  phase  of  the  same  signal  at  some  arbitrarily  chosen  reference 
point.  The  DOA  vector  is  also  referred  to  as  either  the  steering  vector  or  the  array  manifold.  In  the 
case  of  an  array  of  identical  antennas,  the  DOA  vector  depends  only  on  the  azimuth  and  elevation 
angle  of  the  GPS  satellite  and  the  locational  coordinates  of  the  antennas.  In  practice,  though,  the 
antennas  are  not  identical.  In  fact,  due  to  the  limited  aperture  available,  an  array  of  diversely  polarized 
antennas  is  advantageous.  In  addition,  the  limited  aperture  also  necessitates  relatively  close  spacing 
of  the  antennas  thereby  causing  substantial  mutual  coupling  effects.  As  a  result,  the  DOA  vector 
for  a  given  GPS  satellite  is  generally  a  complicated  function  of  many  parameters.  Thus,  a  procedure 
for  estimating  the  DOA  vector  for  a  given  GPS  satellite  is  presented  in  Section  3.4  based  on  the 
cyclostationarity  of  the  C/ A  code  signal. 

How  this  DOA  vector  estimate  is  used  varies  amongst  the  three  advanced  AJ  processors.  It  is 
initially  used  in  beamforming  to  boost  the  Signal-to- White-Noise  Ratio  (SWNR)  for  each  received 
GPS  satellite  signal  by  a  factor  equal  to  the  number  of  antennas,  after  much  of  the  interference  has 
been  canceled  through  power  minimization. 

% 

3.1  Composite  AJ  Processor:  Modular  Stucture  for  Hierarchical  A J/Interference 
Cancellation 

An  embodiment  of  the  fourth  AJ  processor  has  been  developed  that  incorporates  the  first  three  AJ 
processors;  it  entails  three  layers  of  processing  in  a  modular  structure.  Jammers/interferences  are 
canceled  in  a  hierarchical  fashion.  The  first  layer  of  processing  entails  multiple  Space-Time  Power 
Minimization  (ST-PM)  preprocessors  as  presented  in  Section  3.6.  This  first  module  works  to  can¬ 
cel  jammers/interferences  with  respective  power  levels  above  the  noise  floor.  The  second  layer  of 
processing  is  the  per  satellite  combining  of  these  multiple  Space-Time  Power  Minimization  (ST-PM) 
preprocessor  outputs  based  on  an  estimate  of  the  DOA  vector.  This  second  module  works  to  cancel 
jammers/interferers  with  respective  power  levels  near  the  noise  floor.  The  third  layer  of  processing 
entails  a  Minimum  Mean  Square  Error  (MMSE)  receiver,  based  on  the  known  C/A  Gold  Code,  applied 
to  the  output  of  module  2.  This  third  module  works  to  cancel  jammers/interferers  with  respective 
power  levels  below  the  noise  floor,  including  delay  spoofers  and  multipath  propagation.  The  ability  to 
defeat  electronic  attack  in  the  form  of  smart  or  spoofing  jammer  is  a  key  feature  of  the  composite  AJ 
processor. 
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3.2  Key  Enabling  Technology:  Multi-Stage  Nested  Wiener  Filter 

All  of  the  AJ  processors  above  make  use  of  the  Multi-Stage  Nested  Wiener  Filter  (MSNWF)  developed 
by  Dr.  Scott  Goldstein.  The  MSNWF  represents  a  pioneering  breakthrough  in  adaptive  filtering.  It 
provides  rapid  convergence  relative  to  to  standard  LMS  based  adaptations,  and  a  dramatic  reduction 
in  computational  complexity  relative  to'  standard  RLS  based  adaptations,  especially  in  the  case  of 
large  dimension  adaptive  weight  vectors  as  required  for  this  project.  The  MSNWF  is  a  key  enabling 
“technology”  facilitating  real-time  implementation  of  the  space-time  adaptive  algorithms  for  GPS  AJ 
protection  developed  as  part  of  this  project.  A  brief  overview  of  the  MSNWF  is  presented  here. 

Adaptive  filtering  schemes  center  upon  a  linear  Minimum  Mean  Square  Error  (MMSE)  estimation 
problem.  In  any  linear  MMSE  problem,  the  optimum  weight  vector  w  is  the  solution  to  the  Wiener- 
Hopf  equation 

RzzW  =  rdx  (1) 

where  Rxx  is  the  correlation  matrix  of  the  data  and  r^x  is  the  cross-correlation  vector  between  the  data 
and  the  “desired”  signal.  The  LMS  algorithm  or  RLS  algorithm  may  be  used  to  adaptively  compute 
the  solution  to  (1),  allowing  for  the  fact  that  the  second  order  statistics  actually  vary  with  time.  LMS 
is  a  stochastic  gradient  descent  algorithm. 

The  RLS  algorithm  provides  significantly  faster  convergence  than  the  LMS  algorithm  but  at  a  high 
computational  cost,  especially  in  the  case  of  a  “long”  weight  vector.  Adaptation  in  a  high-dimensional 
space  implies  slow  convergence  for  both  the  LMS  and  RLS  algorithms.  Faster  convergence,  can  be 
achieved  by  constraining  the  weight  vector  to  lie  in  a.  low-dimensional  subspace.  The  trick  is  do  this 
with  little  or  no  attendant  increase  in  the  (asymptotic)  mean  square  error  (MSE). 

Prior  work  along  the  lines  of  dimensionality  reduction  restricted  the  weight  vector  w  to  lie  in  a 
subspace  spanned  by  the  Principal  Components  (PC)  or  dominant  eigenvectors  of  RXI.  Although  this 
speeds  up  convergence,  there  is  the  intense  computational  burden  of  computing  the  dominant  eigenvec¬ 
tors  ofRxx.  A  breakthrough  occurred  in  1997  when  Dr.  Goldstein  formulated  the  Multi-Stage  Nested 
Wiener  Filter  (MSNWF);  the  MSNWF  was  first  presented  at  the  IEEE  Military  Communications 
Conference  in  November  1997. 

The  MSNWF  represents  a  pioneering  breakthrough  in  that  it  simultaneously  achieves  a  convergence 
speed-up  substantially  better  than  that  achieved  with  PC  and  a  dramatically  reduced  computational 
burden  relative  to  PC  as  well.  Intuitively  speaking,  achieving  the  best  of  both  worlds  -  faster  con¬ 
vergence  AND  reduced  computation  -  is  made  possible  by  making  use  of  the  information  inherently 
contained  in  both  R^x  and  r^x  in  choosing  the  reduced-dimension  subspace  that  w  is  constrained  to 
lie  within.  In  contrast,  PC  only  makes  use  of  the  information  embedded  in  Rxx. 

In  MSNWF,  there  is  no  computation  of  eigenvectors.  It  has  been  shown  that  MSNWF  implicitly 
constrains  the  weight  vector  to  lie  in  the  Krylov  subspace  spanned  by  {  r^,  R^r^,  R^r*,,  ...  , 
R^r1  rdx  }.  The  word  “implicit”  is  used  since  there  are  implementations  of  the  MSNWF  that  do  not 
require  the  formation  of  the  correlation  matrix.  This  point  is  addressed  later. 

One  key,  although  subtle,  benefit  of  the  MSNWF  is  that  this  algorithm  can  work  in  the  critical 
low-sample  support  operational  environment  where  other  adaptive  algorithms  fail.  In  other  words, 
the  ability  for  rapid  adaptation  is  matched  by  a  lower  requirement  for  training  data  to  estimate  the 
statistics.  Thus  while  many  least-squares  algorithms  or  orthogonal  filter  structures  may  offer  faster 
convergence  than  stochastic  gradient  algorthms,  the  MSNWF  is  the  only  algorithm  which  actually 
reduces  the  sample  support  requirements  without  degrading  performance.  Moreover,  this  presents 
the  ability  to  track  nonstationary  signal  environments.  Finally,  it  has  been  demonstrated  that  the 
MSNWF  is  extremely  robust  to  real-world  sensor  effects  that  lead  to  catastrophic  failure  in  other 
adaptive  algorithms. 

The  MSNWF  algorithm  is  summarized  below.  The  interpretation  of  the  “desired”  signal  do[n] 
varies  amongst  the  different  AJ  processors. 

•  Initialization:  do(n)  and  Xo(n)  =  x(n) 
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Forward  Recursion:  For  k  =  1,2, 


pfc  =.  s{4-i(n)xfc-i(n)}/!i£;{4-i(n)xfc-i(n)}ll 
40)  =  pfxfc-iO) 

B  =  I  -  PkPk 
xfc(n)  =  Bxfc_x(n) 

(2) 


•  Backward  Recursion:  For  k  =  D,D  —  1, 1,  with  ep(n)  =  dp(n): 


=  £J{4_1(n)efc(n)}/JB{|efc(n)|2} 

e/c-iO)  =  4-1 0)  -  WfeCfcO) 


(3) 


It  can  be  easily  shown  that 

—  (I  ~  PfcPf)Rfc-iPfc 

.  P*+i~II(I-ph>?)b*-iP*II 

where 

Rfc+1  =  (I  -  Pfc+iPiH-i)R.fc(I  -  Pfc+iPfe+i)  (5) 

for  k  =  0, 1,  ...,D  —  1,  where  pi  =  r<fx  and  Ro  =  RZi-  It  follows  that  the  matrix  Tp  =  [px  p2  ...  ppj 
contains  orthonormal  columns  and  that  the  reduced  dimension  D  x  D  correlation  matrix  T^RxxTp 
is  tri-diagonal. 

The  MSNWF  is  depicted  in  Figure  1  which  clearly  displays  the  multiple  stages  and  nested  structure. 
Operating  in  a  D-dimensional  space  is  tantamount  to  “cutting  off’  all  stages  below  the  D- th  stage. 
The  updating  of  the  scalar  weights  u>k  in  Figure  1  may  be  effected  through  a  simple  LMS  algorithm. 
We  have  also  have  developed  recursive  adaptations  for  updating  the  so-called  blocking  matrices  Bfc 
in  an  RLS  type  fashion.  The  appropriate  implementation  of  the  MSNWF  for  each  of  the  different 
AJ  processors  is  described  in  the  ensuing  sections.  The  impressive  capabilities  of  the  MSNWF  are 
demonstrated  in  Section  3.3.2  for  the  case  of  a  space-time  power  minimization  based  preprocessor. 


Figure  1:  Structure  of  successive  stages  of  the  Multistage  Nested  Wiener  Filter. 
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3.3  Power  Minimization  Based  Space-Time  Preprocessor 

As  discussed  previously,  the  power  minimization  principle  is.  based  on  the  fact  that  prior  to  correlation, 
the  C/A  code  signal  from  any  GPS  satellite  is  on  the  order  of  30  dB  below  the  noise  floor.  (The  P  code 
is  even  further  below  the  noise  floor.)  A  proven  and  widely  used  means  of  interference  suppression  for 
GPS  is  based  on  the  principle  of  power  minimization.  In  the  case  of  space-only  processing,  the  weight 
for  the  reference  antenna  is  set  to  unity  and  the  weights  for  the1  auxiliary  antennas  are  determined  as 
those  which  drive  the  output  power  of  the  beamformer  as  close  to  the  noise  floor  as  possible.  This 
adaptively  steers  spatial  nulls  in  the  directions  of  those  interferes  whose  power  levels  are  above  the 
noise  floor. 

A  system  based  on  space-only  processing  requires  the  placing  of  a  spatial  null  in  the  direction  of  each 
narrowband  interferer  as  well  as  each  wideband  interferer.  As  discussed  in  Section  2.1,  this  consumes 
precious  degrees  of  freedom  and  potentially  blots  out  angular  sectors  of  space  where  GPS  satellites  lie. 
Again,  both  CRPA  AE  and  GAS- IN  systems  only  allow  for  space- only  processing.  In  contrast,  space- 
time  processing  only  requires  the  formation  of  a  “point-like”  null  in  the  multi-dimensional  spectrum  at 
the  frequency-angle  coordinates  of  each  strong  narrowband  interferer.  A  sharp  line  null,  i.e.,  a  spatial 
null  across  the  entire  frequency  band,  is  still  required  along  the  angular  coordinates  of  each  broadband 
interferer.  These  nulling  characteristics  of  power  minimization  based  space-time  pre-processing  are 
illustrated  in  Figures  3(c)  and  3(f),  which  is  discussed  shortly. 

Since  it  is  desirable  to  track  the  respective  signals  from  as  many  GPS  satellites  as  possible,  the 
space-time  preprocessor  should  “pass”  unaltered  as  many  •  GPS  signals  as  possible  while  canceling 
the  interference.  Mathematically,  it  is  desired  that  the  multidimensional  Fourier  Transform  of  the 
space-time  weights  be  as  flat  in  magnitude  as  possible.  In  the  case  of  space-time  power  minimization, 
distortion  to  a  GPS  signal  from  a  given  satellite  can'  occur  if  a  narrowband  interferer  arrives  from 
the  same  azimith-elevation  direction  as  that  satellite.  Note,  though,  that  if  one  employs  space- 
only  processing  to  take  out  a  narrowband  interferer  under  these  same  conditions,  the  particular  GPS 
satellite  is  taken  out  by  the  spatial  null  as  well. 

The  disadvantage  of  space- time  processing  relative  to  space-only  processing  is  the  large  dimension¬ 
ality  of  the  space-time  weight  vector.  Referring  to  Figure  2,  there  are  N  taps  per  each  of  M  antennas 
ostensibly  requiring  adaptation  in  an  IVAf-dimensional  space,  as  opposed  to  an  Ilf-dimensional  space 
in  the  case  of  space-only  processing.  The  higher  dimensionality  translates  into  a  large  computational 
burden  and  slow  convergence.  The  MSNWF  is  employed  to  counter  these  deleterious  effects.  Il¬ 
lustrative  preliminary  simulation  results  are  presented  in  Section  3.3.2  demonstrating  the  impressive 
capabilities  of  the  MSNWF. 

The  primary  innovation  in  the  space-time  extension  of  the  power  minimization  based  preprocessor, 
developed  as  part  of  this  effort  as  a  basic  module  for  the  more  advanced  AJ  processors  presented,  is 
effecting  it  adaptively  with  the  MSNWF.  There  are  additional  advantages  over  and  above  the  speed¬ 
up  in  convergence  and  reduction  in  computational  burden  afforded  by  the  MSNWF.  For  example,  it 
allows  for  a  larger  number  of  “taps”  per  antenna  (large  N )  without  taking  too  much  of  a  “hit”  in 
terms  of  convergence  and  computational  burden.  This  facilitates  the  formation  of  sharper  spectral 
nulls  in  canceling  narrowband  interferers,  implying  less  distortion  to  the  GPS  signals. 

3.3.1  Formulation  of  Objective  Function 


The  data  is  obtained  by  baseband  I  &  Q  frequency  conversion  at  each  antenna  followed  by  A/D 
conversion  at  a  Nyquist  rate  commensurate  with  the  20.46  MHz  band  occupied  by  the  P  code  signal. 
Referring  to  Figure  2,  xm(n)  is  an  N  x  1  vector  containing  N  successive  samples  of  the  output  of  the 
m-th  antenna  sampled  at  a  rate  above  or  equal  to  the  Nyquist  rate  for  the  P(Y)  signal. 

xm(n)  =  [xm(n),xm(n-  l),.:.,xm(n-  N  +  1)}T,  m  =  l,2,...,M.  (6) 
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Figure  2:  Adaptive  filter  structure  for  power  minimization  based  space-time  preprocessor. 


Similarly,  the  N  tap  weights  for  the  m-th  antenna  are  placed  as  the  components  of  an  N  x  1  vector  as 

wm  =  [wm(0),wm(l),...,wm(N  -  1)]T,  m  =  1,2, M.  (7) 

The  NM  x  1  space-time  snapshot,  x(n),  and  the  NM  x  1  weight  vector,  w  are  formed  from  concate¬ 
nating  xm(n),  m  =  1,2, ...,  M,  and  wm,  m  =  1, ...,  M,  respectively,  as 

xi(n)  1  W! 

x2(n)  ,  w2 

w  = 

Xm(")  J  wM 

To  incorporate  the  unity  weight  constraint  on  the  first  tap  of  the  m-th  antenna,  we  define  x_m(n) 
as  the  ( NM  —  1)  x  1  sub- vector  of  x(n)  containing  all  but  the  first  element  of  xm(n),  xrn(n).  Similarly, 
\v  is  defined  as  the  (NM  —  1)  x  1  sub-vector  of  w  containing  all  but  the  first  element  of  wm.  With 
these  definitions,  the  power  at  the  preprocessor  output  may  be  expressed  as 

£  {|wHx(n)|2}  =  £{\  xm(n)  +  wHx_m(n)|2}.  (9) 

where  £  denotes  expected  value  in  a  statistical  sense  and  (■')n  denotes  conjugate  transpose.  Expressing 
the  preprocessor  output  power  in  this  fashion  facilitates  an  adaptive  filtering  formulation  where  the 
output  of  the  first  tap  of  the  m-th  antenna  serves  as  the  “desired”  signal  such  that  the  “error”  signal  is 
xm(n)  +  wFx_m(n)  .=  xm(ri)  —  (—wH^m(n)).  This  facilitates  the  use  of  LMS  and/or  RLS  based  per 
sample  adaptations  in  which  the  “desired  signal”  is  the  output  of  the  first  tap  at  the  m-th  antenna. 
For  the  MSNWF,  do(n)  =  xm(n)  in  Figure  1. 
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Jammer  Type 

Power 

AO  A 
Ex.l 

AOA 
Ex. 2 

Bandwidth 

Wideband 

-100  dBW 

20° 

20° 

20  MHz 

Wideband 

-110  dBW 

0° 

20  MHz  . 

Wideband 

-100  dBW 

— 

1 

to 

o 

o 

20  MHz 

Wideband 

-100  dBW 

— 

o 

O 

1 

20  MHz 

Wideband 

-110  dBW 

-60° 

20  MHz 

Jammer  Type 

SNR 

AO  A 

AOA 

Freq.  relative  to  LI 

Narrowband 

-100  dBW 

60° 

o 

O 

SO 

-10  MHz 

Narrowband 

-100  dBW 

15° 

— 

-5  MHz 

Narrowband 

-100  dBW 

-10° 

— 

0  MHz 

Narrowband 

-100  dBW 

o 

O 

CO 

1 

- 

5  MHz 

Narrowband 

-110  dBW 

-55° 

_ 

10  MHz 

Table  1:  Interference  parameters  for  illustrative  simulations  presented  in  Figure  3.  Noise  floor  is  at 
-127  dBW. 

3.3.2  Illustrative  Proof-of-Concept  Simulation  Examples:  MSNWF  Implementation  of 
Power  Minimization  Space-Time  Preprocessor 

The  illustrative  proof-of-concept  simulation  result  presented  here  employed  N  =  7  taps  at  each  of 
M  —  7  antenna  elements  equi-spaced  along  a  line.  Although  typical  antenna  arrays  for  GPS  are  two- 
dimensional  (planar),  circular,  for  example,  or  conformal,  a  linear  array  was  used  in  this  illustrative 
simulation  example  in  order  to  have  only  one  angular  variable  for  display  purposes.  This  allows  the  use 
of  a  single  mesh  or  contour  plot  to  display  the  two-dimensional  Fourier  Transform  of  the  space-time 
weights  obtained  from  a  given  run. 

In  addition,  N  =  7  is  a  very  small  number  of  taps  to  employ  at  each  antenna.  A  significantly  larger 
number  is  needed  in  practice  in  order  to  form  sharp  ‘point-nulls’  at  the  angle-frequency  coordinates 
of  a  narrowband  interferer  and  thereby  minimize  distortion  to  the  GPS  signal.  Finally,  the  number  of 
jammers  employed  in  the  simulations  was  arbitrarily  chosen  to  be  six.  The  number  of  jammers  that 
can  be  handled  in  practice  depends  on  the  type  of  jammers  and  varies  amongst  the  various  AJ  algo¬ 
rithms  with  the  more  complex  algorithms  being  able  to  handle  a  larger  number  of  interferers.  Again, 
the  results  presented  here  are  merely  preliminary  proof-of-concept  illustrative  simulations  providing 
evidence  of  the  power  of  the  MSNWF  algorithm. 

Table  1  summarizes  the  parameters  of  the  wideband  and  narrowband  jammers  simulated  in  each' 
of  two  different  scenarios.  The  angles-of-arrival  (AOA’s)  listed  are  relative  to  broadside;  the  frequency 
listed  for  each  narrowband  jammer  is  its  offset  relative  to  LI.  The  first  scenario  involved  five  narrow- 
band  jammers  and  one  wideband  jammer.  The  attendant  results  are  plotted  in  Figure  3(a)  through 
3(c).  The  second  scenario  involved  one  narrowband  jammer  and  five  wideband  jammers.  The  atten¬ 
dant  results  are  plotted  in  Figure  3(d)  through  3(f).  In  all  cases,  the  first  tap  at  the  first  antenna  was 
constrained  to  be  unity.  Given  the  20.46  MHz  receiver  bandwidth  at  each  antenna,  the  noise  floor 
was  determined  to  be  approximately  -128  dBW.  Recall  the  goal  of  power  minimization  is  to  drive  the 
output  power  of  the  space-time  beamformer  as  close  to  the  noise  floor  as  possible. 

Convergence  Performance.  Figures  3(a)  and  3(d)  plot  the  average  power  output  of  the  space- 
time  power  minimization  preprocessor  based  on  the  MSNWF  as  .  a  function  of  subspace  dimension 
or  rank  of  the  dimensionality  reducing  matrix  transformation.  The  subspace  dimension  at  which 
MSNWF  approximately  achieves  the  performance  of  the  full-dimension  ideal  (asymptotic)  Wiener 
filter  is  roughly  the  same  in  both  scenarios,  around  eight.  In  contrast,  Principal  Components  (PC) 


Baseband  Frequency  (MHz)  w  Power  (dBW)  Power  (dBW) 


Rank  Rank 


3(a).  Power  Output  vs  Rank:  1  WB,  5  NB.  3(d).  Power  Output  vs  Rank:  5  WB,  1  NB. 


Average  Power  Output  of  Preprocessor  Average  Power  Output  of  Preprocessor 


Space-time  Snapshots  Space-time  Snapshots 


..  Power  Output  vs  Snapshots:  1  WB,  5  NB.  3(e).  Power  Output  vs  Snapshots:  5  WB,  1  NB. 


Theta  (Angle  of  Arrival)  Theta  (Angle  of  Arrival) 

3(c).  Angle-Frequency  Response:  1  WB,  5  NB.  3(f).  Angle-Frequency  Response:  5  WB,  1  NB. 
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generally  requires  a  subspace  dimension  equal  to  the  number  of  degrees  of  freedom  taken  up  by  the 
jammers  to  achieve  the  same  output  power  level.  Each  narrowband  jammer  takes  up  one  degree  of 
freedom.  Each  wideband  jammer  takes  up  N  —  7  degrees  of  freedom,  where  N  is  the  number  of 
taps  per  antenna.  This  is  because  the  cancellation  of  a  wideband  jammer  requires  a  spatial  null, 
implying  a  null  across  the  entire  20.46  MHz  spectrum  at  its  AOA.  In  Scenario  1,  the  jammers  take  up 
5x1  +  1x7=12  degrees  of  freedom;  in  Scenario  2,  the  jammers  take  uplxl  +  5x7  =  36  degrees 
of  freedom.  1 

Figures  3(b)  and  3(e)  examine  the  space-time  snapshot  sample  support  necessary  to  effectively  null 
the  jammers  for  each  of  the  two  scenarios  simulated.  The  power  output  for  each  sample,  support  level 
was  averaged  over  250  Monte  Carlo  trial  runs.  The  greatest  differential  in  performance  between  the 
MSNWF  and  PC  based  methods  is  observed  in  Figure  3(e)  corresponding  to  Scenario  2.  In  this  case, 
Figure  3(d)  and  the  above  calculation  dictate  that  PC  needs  to  adapt  in  a  36-dimensional  subspace, 
while  the  MSNWF  need  only  adapt  in  a  10-dimensional  space.  This  allows  MSNWF  to  converge  more 
rapidly  than  PC. 

Computational  Performance.  At  the  same  time,  there  is  no  computation  of  eigenvectors  involved 
in  the  MSNWF.  Only  a.  small  number  of  simple  matrix- vector  multiplications  are  needed  to  determine 
the  required  low-dimensional  subspace.  In  contrast,  the  computation  of  the  36  eigenvectors  needed 
by  the  PC  method  in  Scenario  2  is  a  substantial  computational  burden.  Note,  though,  that  PC  does 
indeed  converge  more  quickly  than  the  full-dimension  (49),  finite-sample  Wiener  filter.  Note  that 
Figure  3  also  displays  the  performance  of  the  Cross-Spectral  Metric  (CSM)  method.  Similar  to  the 
PC  method,  the  CSM  method  constrains  the  space-time  weight  vector  to  lie  in  a  subspace  spanned  by 
a  subset  of  the  eigenvectors  of  the  space-time  correlation  matrix.  The  choice  of  eigenvectors  is  dictated 
by  a  cross-spectral  metric  derived  from  the  cross-correlation  vector,  rather  than  simply  choosing  those 
eigenvectors  associated  with  the  largest  eigenvalues.  CSM  yields  improved  performance  relative  to . 
PC,  but  its  performance  is  not  nearly  as  good  as  MSNWF  and  it  too  requires  the  computation  of 
eigenvectors. 

Nulling  Performance/Distortion  Issues.  Figures  3(c)  and  3(f)  display  contour  plots  of  the  magni¬ 
tude  of  the  multi-dimensional  Fourier  Transform  of  the  space-time  weights  obtained  from  the  MSNWF 
with  40  space-time  snapshots.  For  Scenario  1,  Figure  3(c)  displays  a  well-defined  “point-null”  at  the 
angle- frequency  coordinate  of  each  narrowband  jammer  and  a  well-defined  “line-null”  along  the  arrival 
angle  of  the  wideband  jammer.  For  Scenario  2,  Figure  3(f)  displays  a  well-defined  “point-null”  at  the 
angle- frequency  coordinate  of  the  one  narrowband  jammer  and  a  well-defined  “line-null”  along  the 
respective  arrival  angle  of  each  of  the  five  the  wideband  jammers.  As  important,  in  both  cases  the 
response  of  the  space-time  beamformer  is  observed  to  be  relatively  flat  away  from  the  null  locations. 

3.4  DOA  Vector  Estimation  Based  on  Cyclostationarity  of  C/A  Code  Signal 

The  other  three  AJ  processors  developed  as  part  of  this  effort  produce  a  different  space-time  filtered 
output  for  each  GPS  satellite  in  the  FOV.  A  key  component  of  each  of  these  three  AJ  processors 
is  exploitation  of  the  cyclostationarity  of  the  C/ A  code  signal  to  estimate  the  Direction-of- Arrival 
(DOA)  vector  for  a  given  GPS  satellite.  Again,  the  m-th  element  of  the  M  x  1  DOA  vector  for  a 
given  GPS  satellite  is  the  gain  and  phase  of  the  GPS  signal  received  from  that  satellite  at  the  m-th 
antenna,  relative  to  the  gain  and  phase  of  the  same  signal  at  some  arbitrarily  chosen  reference  point. 
Exploitation  of  cyclostationarity  enables  one  to  estimate  the  DOA  vector  for  a  GPS  satellite  even  in 
the  presence  of  strong  interference  from  multiple  sources. 

How  this  DOA  vector  estimate  is  used  varies  amongst  the  different  methods.  It  is  initially  used  in 
beamforming  to  boost  the  Signal-to- White-Noise  Ratio  (SWNR)  for  each  received  GPS  satellite  signal  • 
by  a  factor  equal  to  the  number  of  antennas,  after  much  of  the  interference  has  been  canceled  through 
power  minimization.  This  is  discussed  shortly. 

The  DOA  vector  obtained  by  processing  the  C/A  code  signal  may  be  used  advantageously  in 
processing  the  P  code  signal  from  that  same  satellite  since  the  two  signals  arrive  from  the  same 
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Figure  5:  Cyclostationarity  based  method  for  estimating  the  DOA  vector  for  the  k-th  GPS  satellite 
in  the  field  of  view.  As  a  result  of  Doppler  effects,  the  block,  shown  within  a  dashed  box  (and  each 
of  the  other  same  blocks)  is  implemented  via  a  bank  of  parallel  filters  effecting  partial  despreading  as 
shown  in  in  Figure  6. 


(n  •  (A 


Figure  6:  Partial  despreading  via  a  bank  of  L  parallel  filters  derived  as  length  Q  (nonoverlapping) 
sub-blocks  of  the  1023  length  Gold  Code  for  the  k-th  GPS  satellite. 


direction.'  This  is  discussed  later.  The  C/A  code  signal  is  used  to  find  the  DOA  vector  for  a  given 
GPS  satellite  since  it  is  cyclostationary,  arising  from  the  fact  that  the  same  1023  length  Gold  Code 
is  repeated  every  millisecond.  The  P  code  signal  is  cyclostationary  but  on  a  much,  much  longer  time 
scale. 

3.4.1  Partial  Despreading  Due  to  Doppler  Shift  Effects 

Estimation  of  the  DOA  vector  for  a  given  GPS  satellite  is  facilitated  by  boosting  its  SWNR  and  its  * 
SINR  as  much  as  possible,  prior  to  cyclic  spatial  correlation  matrix  analysis  (to  be  explained  shortly). 
To  this  end,  we  first  perform  digital  lowpass  filtering  at  each  antenna,  as  indicated  in  Figure  5,  since  the 
C/A  code  signal  only  occupies  one-tenth  of  the  bandwidth  occupied  by  the  P  code  signal.  Accordingly, 
the  filtered  signal  at  each  antenna  may  be  decimated  or  downsampled  by  a  factor  of  ten  to  produce 
the  outputs,  ym(n),  m  =  1, ...,  M.  The  combination  of  the  lowpass  filtering  and  decimation  has  been 
efficiently  implemented  through  multirate  DSP  using  polyphase  filters. 

Ideally,  this  would  be  followed  by  correlating  and  summing  with  the  C/A  code  for  the  fc-th  GPS 
satellite  in  the  FOV,  sliding  over  one  chip  (over  one  half-chip  with  two  times  per  chip  sampling  of  C/A 
signal)  at  a  time.  As  indicated  in  Figure  5,  continual  correlation  in  this  fashion  is  equivalent  to  running 
each  output  ym(n)  through  a  filter  whose  impulse  response  is  the  time-reverse  of  the  C/A  code.  Since 
the  C/A  code  is  repeated  twenty  times  in  one  bit  of  navigation  information,  an  additional  13  dB  boost 
in  SWNR  could  be  achieved  by  continuous  correlation  with  a  sequence  formed  from  twenty  repetitions 
of  the  C/A  code  for  the  fc-th  satellite.  However,  the  signal  arriving  from  each  GPS  satellite  may  be 
Doppler  shifted  as  much  as  up  to  dr  5  KHz  away  from  LI  (or  L2)  due  to  motion  of  the  satellite,  thereby, 
effectively  destroying  the  processing  gain. 

In  fact,  the  Doppler  shifts  are  potentially  large  enough  that  the  variation  in  phase  over  just  one 
Gold  Code  block  is  substantial.  Now,  the  cyclostationarity  based  method  for  estimating  the  DOA 
vector  is  immune  to  Doppler  effects  as  a  result  of  both  an  inherent  magnitude  squaring  and  the  fact 
that  the  Doppler  shift  for  a  given  GPS  satellite  is  the  same  at  each  receive  antenna.  However,  its 
efficacy  depends  on  the  SWNR  of  the  satellite  signal  of  interest  being  “boosted”  above  other  signals 
that  repeat  at  the  same  rate. 

This  is  where  the  partial  despreading  is  used.  The  idea  is  to  break  the  length  1023  Gold  Code 
for  a  given  GPS  satellite  into  L  nonoverlapping  sub-blocks  of  length  Q.  The  l- th  sub-block  of  the  k- th 
Gold  Code  c^k\n)  is  denoted  c^k\n\  t).  Q  is  chosen  as  large  as  possible  under  the  constraint  that  the 
change  in  phase  over  the  sub-block  due  to  Doppler  be  small  enough  to  achieve  a  processing  gain  for 
each  sub-block  as  close  to  Q  as  possible.  Since  interfering  signals  that  do  not  repeat  at  the  same  rate 
as  the  C/A  code  asymptotically  have  vanishing  contribution  to  the  cyclic  spatial  correlation  matrix  (to 
be  mathematically  defined  shortly),  a  processing  gain  of  18  dB  is  more  than  sufficient  for  the  purpose 
of  estimating  the  DOA  vector  for  the  fc-th  satellite  via  the  cyclostationarity  scheme  presented  in  the 
next  section. 

In  the  trade-off  between  processing  gain  and  Doppler  effects,  a  reasonable  possibility  for  starters 
is  to  decompose  the  Gold  Code  into  L  =  16  blocks  of  length  Q  =  64.  For  the  sake  of  notational 
simplicity,  the  last  block  is  padded  with  a  single  zero.  Optimal  designs  have  been  developed  though 
performance  analysis. 

Thus,  as  depicted  in  Figure  6,  the  outputs  ym(n)  are  then  passed  through  a  bank  of  L  parallel 
filters  where  the  impulse  response  of  the  f-th  filter  is  the  time-reverse  of  the  f-th  sub-block  of  the  Gold 
Code  associated  with  the  fc-th  GPS  satellite,  denoted  c^k\n\t).  This  creates  the  outputs  Zm\n\£), 
where  m  is  the  antenna  index,  k  is  the  satellite  index  and  i  is  the  sub-block  index. 

3.4.2  Cyclic  Correlation  Matrix  Analysis 

The  above  process  yields  a  “finger”  in  .the  Mh  sub-block  output  at  the  m-th  antenna  every  1023 
samples,  multiplied  by  both  a  navigation  information  bit  value  and  the  m-th  value  of  the  DOA  vector. 
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At  this  point,  the  relative  locations  of  the  fingers  are  unknown  and  they  may  not  be  detectable  due 
to  interference  despite  the  18  dB  processing  gain.  For  the  £-th  sub-block,  we  thus  propose  to  form  the 
following  cyclic  spatial. correlation  matrix  of  dimension  equal  to  the  number  of  antennas,  M. 

R *z(£)  e“j27r k?23  z(fc) (n;  (n;  £)  (10) 

Ns  n= o 

where  z W(n)  =  [z[k\n;£),z^k\n;e)z$ {n\l))T  and  Ns  is  the  number  of  C/A  code  blocks  over  which 
R^(^)  is  averaged.  Note  that  the  time  duration  over  which  Ra  is  averaged  may  span  numerous 
cycles  of  the  Doppler  shift  modulation  since  the  asymptotic  structure  of  Ra  is  not  affected  by  Doppler. 
The  immunity  of  Ra  to  Doppler  effects  is  a  result  of  both  the  magnitude  squaring  inherent  in  the 
conjugate  transpose  operation  denoted  by  H  and  the  fact  that  the  Doppler  shift  is  the  same  at  each 
receive  antenna. 

In  the  asymptotic  limit,  IVfz(£)  only  has  contributions  from  signals  that  have  the  same  cycle 
frequency,  repeat  every  1023  samples.  Any  jammer  that  does  not  repeat  at  this  rate  will  asymptotically 
have  vanishing  contribution  to  R %z(£),  as  does  the  noise.  In  the  case  of  no  delay  spoofers,  we  have 

shown  that  the  dominant  contribution  to  Ra(^)  is  proportional  to  e~^2lx  1023“ where  a®  is 
the  DOA  vector  for  the  k- th  GPS  satellite  and  rk  is  the  relative  time  location  of  the  “finger”  in  the 
first  sub-block  code  filter  output,  relative  to  a  single  period  of  the  C/A  code  spanning  1023  chips. 

The  phase  factor  accounts  for  the  fact  that  the  “fingers”  occurring  in  the  ^-th  sub-block 

output  are  offset  in  time  by  £Q  units  relative  to  the  “fingers’"  occurring  in  the  first  sub-block  output, 

Z  =  0, 1, ...,  L  —  1. 

Two  methods  for  averaging  the  L  values  of  have  been  investigated.  First,  they  may  be 

incoherently  averaged  as 

<n> 

^=0 

Recall  that  the  dominant  contribution  to  Ra(^)  is  proportional  to  e_j27r  1023  which  is  not 

Hermitian-symmetric.  Second,  the  L  values  of  R“z(£)  may  be  coherently  averaged  as 

=  (12) 
L  £— 0 

The  weighting  A  is  ideally  of  the  form  e~j  A ,  where  A  is  the  center  phase  offset  of  one  sub-block  relative  ■ 
to  the  next  due  to  the  Doppler  modulation.  Procedures  for  estimating  A  to  elevate  the  contribution  of 
the  desired  satellite  signal,  over  and  above  the  18  dB  processing  gain  achieved  with  partial  despreading 
with  sub-blocks  of  length  Q  =  64,  have  been  developed  as  part  of  this  project. 

Note  that  Section  3.6  proposes  applying  this  DOA  estimation  procedure  at  the  output  of  a  bank  of 
parallel  ST- PM  preprocessors,  as  shown  in  Figure  9,  where  interferers  above  the  noise  floor  are  already 
canceled.  In  this  case,  only  jammers  below  the  noise  floor  that  repeat  at  this  rate  will  potentially 
contribute  significantly  to  RJZ.  Given  that  the  contribution  from  the  A;-th  GPS  satellite  benefits  from 
an  18  dB  processing  gain,  the  only  types  of  jammers  that  will  contribute  significantly  to  RJZ  are  delay 
spoofers. 

In  the  case  of  no  delay  spoofers,  the  dominant  contribution  to  Ha  is  proportional  to 
Thus,  a  couple  of  power  iterations  may  be  employed  to  estimate  the  dominant  left  singular  vector  of 
Ra  and  thereby  provide  an  estimate  of  the  DOA  vector  aSk )  (to  within  an  unknown  multiplicative 
scalar  inconsequential  for  beamforming  purposes). 

If  there  are  J  delay  spoofers,  R %  will  be  of  the  form 

J 

=  e^27T  1023  aS^aS^*1  +  e-j27rio23  (13) 

p= i 
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where  is  the  DOA  vector  for  the  g-th  delay  spoofer  and  rp  is  its  corresponding  relative  time  delay. 
Methods  for  canceling  delay  spoofers  are  presented  shortly. 

3.4.3  Proof-of-Concept  Experiment:  Cyclostationarity  Based  Interference  Cancellation 

The  efficacy  of  expioiting  the  cyclostationarity  of  the  C/A  code  signal  to  estimate  the  DOA  vector  for 
a  given  GPS  satellite  and  ultimately  combat  interference  is  demonstrated  using  data  derived  from  an 
experimental  test-bed  built  at  the  Polytechnic  University  of  Madrid  through  a  collaboration  between 
Dr.  Zoltowski  and  the  Electromagnetics  Group  there.  A  brief  description  of  the  test-bed  is  provided 
below. 

A  block  diagram  is  shown  in  Figure  7  (a).  The  antenna  array  consists  of  6  elements  equi-spaced 
along  a  circle  of  radius  equal  to  roughly  half- wavelength  at  LI.  Each  antenna  element  is  a  low 
cost  stacked  patch  antenna.  This  configuration  provides  a  wide  bandwidth  allowing  the  reception 
of  INMARSAT  as  well  as  GPS.  Tp  achieve  circular  polarization,  a  90°  hybrid  is  placed  below  each 
element.  (A  picture  of  the  experimental  prototype  in  Figure  8  (a)  shows  the  physical  antenna  array 
structure.)  Following  this  are  six  Low  Noise  Amplifiers  (LNA)  in  MMIC.  The  RF  outputs  of  the  six 
antennas  are  handed  to  six  standard  DS-SS  receivers  with  the  special  feature  that  all  of  them  are  locked 
to  the  same  Local  Oscillator  (LO).  These  six  spread  spectrum  receivers  are  basically  composed  of  two 
stages:  downconversion  from  RF  to  IF,  followed  by  digital  down  conversion  from  IF  to  BaseBand 
(BB).  Only  the  C/A  code  signal  was  processed  in  this  prototype. 

An  experiment  was  conducted  in  which  the  experimental  GPS  receiver  was  programmed  to  receive 
a  GPS  signal  arriving  at  an  elevation  angle  of  30  degrees  with  respect  to  the  boresite  axis  of  the 
array.  After  despreading,  the  SNR  of  the  GPS  signal  was  roughly  10  dB  per  element.  In  addition, 
interference  was  intentionally  injected  by  a  nearby  radiating  antenna  at  an  angle  of  20  degrees  with 
respect  to  boresite,  and  at  a  power  level  40  dB  above  the  desired  GPS  signal  (prior  to  despreading). 
A  6  x  6  cyclic  spatial  correlation  matrix,  Rzz,  was  formed  from  the  outputs  of  each  antenna  in 
the  manner  described  previously.  In  addition,  a  standard  spatial  correlation  matrix  was  formed  as 
=  4  £&'*“>(  n)z^H(n).  The  interferer  did  not  have  the  same  cycle  frequency  as  the  C/A 
code.  Therefore,  Rzz  essentially  had  one  dominant  term  proportional  to  sSk^aSk^H.  The  standard 
correlation  matrix  was  approximately  of  the  form  Hzz  —  +  apa^a.^H  +  aw I,  where  crjfc, 

■<j/,  and  aw  reflect  the  relative  power  levels  of  the  jammer,  the  GPS  signal  (post-correlation),  and  the 
receiver  noise,  respectively,  is  the  DOA  vector  for  the  interferer. 

Rather  than  compute  an  estimate  of  the  GPS  DOA  vector  a ^  from  the  largest  left  singular  vector 
of  R“z  (since  the  array  was  not  well-calibrated),  a  set  of  beamforming  weights  was  determined  by 
computing  the  largest  generalized  eigenvector  of  the  matrix  pencil  {RZZ,R22:}.  In  order  to  give  some 
idea  of  the  improvement  achieved,  in  Figure  7  (c)  the  eye  diagram  at  the  antenna  element  with  the 
best  SNR  is  compared  with  that  obtained  at  the  beamformer  output.  It  is  clear  that  the  jammer  has 
been  sufficiently  rejected  so  as  to  “open  up  the  eye.” 

A  very  important  aspect  of  exploiting  cyclostationarity  is  that  it  steers  the  beam  and  adapts  the 
nulls  based  only  on  the  known  C/A  code.  Due  to  nonlinear  effects  amongst  the  individual  antennas, 
the' most  predominant  of  which  arise  from  mutual  coupling,  it  is  difficult  and  expensive  to  design 
antenna  arrays  calibrated  well  enough  to  achieve  beamforming  through  phase  shifts  derived  from  the 
azimuth  and  elevation  coordinates  of  a  given  GPS  satellite  and  the  locations  of  the  antennas. 

3.5  Per  Satellite  Power  Minimization  Space-Time  Preprocessor 

The  Space-Time  Power  Minimization  (ST-PM)  based  preprocessor  does  NOT  yield  maximal  SINR  for 
any  GPS  satellite,  but  rather  attempts  to  “pass”  all  GPS  satellite  signals  in  the  FOV  as  undistorted 
as  possible  while  canceling  the  interference.  Two  shortcomings  of  the  ST-PM  based  preprocessor  are 
(1)  it  does  not  attempt  to  minimize  distortion  to  any  one  GPS  satellite  signal  and  (2)  it  does  not 
take  advantage  of  the  SWNR  gain  possible  through  co-phasal  beamforming  for  a  given  GPS  satellite. 
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(b)  Structure  of  IF-BB  Stage. 
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(a)  Block  diagram  of  prototype/test-bed.  (c)  Eye  diagrams  from  experimental  data. 


Figure  7:  Experimental  prototype  developed  at  the  Polytechnic  University  of  Madrid  in  collaboration 
with  Purdue  University. 

It  is  presented  that  the  estimate  of  the  DOA  vector  for  a  given  GPS  satellite,  a*,,  obtained  via  the 
presented  cyclostationarity  based  scheme,  be  used  to  achieve  these  goals  as  follows. 

First,  a  different  ST-PM  preprocessor  is  needed  for  each  of  the  K  GPS  satellites  in  the  FOV.  Thus, 
achieving  the  above  two  goals  to'  enhance  GPS  protection  comes  at  a  significant  computational  cost. 
Efficient  implementation  of  the  K  parallel  constrained  ST-PM  preprocessors  using  the  MSNWF  have 
been  developed  as  part  of  this  project.  Since  a  different  space-time  weight  vector  is  formed  for  each 
satellite,  this  is  denoted  by  placing  a  superscript  k  on  the  weight  vector  w. 

For  the  k- th  satellite,  we  constrain  the  inner  product  between  and  the  Mx  1  vector  of  weights  as¬ 
sociated  with  the  same  time  instant  but  spanning  the  M  antennas,  =  [w^  ( n),w ^ (n), ...,  wjffl  (n)]T, 
to  be  unity  for  each  of  the  N  “tap  times”  comprising  the  space-time  adaptive  filter  structure.  This 
leads  to  power  minimization  with  N  linear  constraints: 

Minimize  (14) 

Subject  to:  w^HTnkk  =  1,  n  =  0, 1, ...,  N  —  1 

where  T n  is  a  selection  matrix  composed  of  ones  and  zeroes  such  that  =  TR  w1-^,  where  wjf'  = 
[«)f  (n))ti)^(n))...lii;^(n)]T,  as  defined  above.  Accommodating  the  N  linear  constraints  consumes 
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(a)  Prototype  stacked  patch  antenna  array.  (b)  Compact  array  for  DF  at  HF. 


Figure  8:  Examples  of  experimental  compact  arrays  with  polarization  diversity. 

N  out  of  NM  degrees  of  freedom.  Implementations  of  the  MSNWF  that  accommodate  multiple  linear 
constraints  have  been  developed  as  part  of  this  project. 

Again,  a  different  constrained  ST-PM  filtered  output  is  created  for  each  satellite  in  the  FOV.  If  a 
retro-fit  is  required  to  existing  GPS  UE,  the  outputs  obtained  by  doing  this  for  each  of  the  K  satellites 
in  the  FOV  may  be  simply  summed  and  fed  to  the  GPS  receiver  (after  D/A  and  frequency  conversion 
back  to  LI).  This  is  suboptimal  as  well  as  inefficient.  If  the  AJ  functions  were  integrated  into  the  GPS 
receiver  and  operated  synergistically  with  its  standard  functions,  the  1-D  signal  created  for  the  fc-th 
GPS  satellite  created  by  this  method  would  be  sent  directly  to  a  code-tracker  (and  Doppler  tracker) 
matched  to  that  satellite. 

3.6  Multiple  Power  Minimization  Space-Time  Preprocessors  and  Per  Satellite 
i  Combining 

Although  the  procedure  in  Section  3.4  for  estimating  the  DOA  vector  for  the  fc-th  satellite  is  effective 
even  in  the  presence  of  strong  interference,  it  may  require  averaging  over  many  bits  before  the  interfer¬ 
ence  contribution  to  the  cyclic  correlation  matrix  is  small  enough  to  obtain  a  reliable  estimate  of  a^) . 
We  thus  propose  a  methodology  for  canceling  as  much  interference  as  possible  before  estimating  the 
DOA  vector.  The  architecture  presented  for  such  may  be  used  as  the  first  stage  of  a  multi-stage  AJ 
process,  where  the  first  stage  attempts  to  cancel  jammers  above  the  noise  floor.  This  facilitates  can¬ 
cellation  of  other  AJ  waveforms  hovering  at  or  below  the  noise  floor  in  subsequent  stages  as  discussed 
in  Section  3.7. 

The  key  idea  is  to  create  M  different  ST-PM  preprocessors  as  shown  in  Figure  9:  the  m-th  one  is 
created  by  constraining  the  weight  value  for  the  first  tap  at  the  m-th  antenna  to  be  unity;  the  remaining 
space-time  weights  are  chosen  as  those  which  yield  minimum  output  power.  Efficient  implementations 
of  the  M  parallel  ST-PM  preprocessors  using  the  MSNWF  have  been  developed. 

This  yields  M  linearly  independent  ST-PM  outputs- as  shown  in  Figure  9.  In  each  case,  interferers 
above  the  noise  floor  are  ideally  canceled.  The  magnitude  of  the  angle- frequency  spectrum  of  the  m-th 
set  of  space-time  weights  is  relatively  smooth  with  point  nulls  at  the  angle-frequency  coordinates  of 
each  narrowband  jammer,  and  a  line  null  in  frequency  at  the  angle  of  each  broadband  jammer.  The  idea 
is  to  use  the  estimated  DOA  vector  for  the  k-th  GPS  satellite  to  combine  the  M  ST-PM  outputs  so  that 
the  response  of  the  angle-frequency  spectrum  of  the  composite  sum  of  weighted  ST-PM  preprocessors , 
the  m-th  one  weighted  by  a®  as  shown  in  Figure  9,  is  elevated  along  the  angular  direction  of  that 
satellite,  i.e.,  so  that  its  SWNR  is  ideally  boosted  by  a  factor  of  M.  This  is  achieved  through  co-phasal 
beamforming  where  otrn  is  equal  to  the  complex  conjugate  of  the  DOA  vector  associated  with  the  k-th 
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GPS  satellite. 

The  M  different  ST-PM  weights  transform  the  data  from  the  original  M-dimensional  element-space 
to  an  M-dimensional  reduced  interference  beamspace.  The  DOA  vector  for  each  satellite  will  generally 
be  altered  through  this  transformation  process.  It  is  thus  necessary  to  apply  the  cyclostationarity 
based  DOA  vector  estimation  module  in  Figure  5  at  the  outputs  of  the  M  ST-PM  preprocessors.  In 
fact,  it  is  highly  advantageous  to  do  this  since  the  interference  is  substantially  reduced  at  the  outputs 
of  the  M  ST-PM  preprocessors.  Combined  with  the  processing  gain  achieved  by  correlating  and 
summing  with  sub-blocks  of  the  C/A  code  for  the  k- th  satellite,  the  reduction  in  interference  causes 
the  “largest”  left  singular  vector  of  the  M  x  M  cyclic  correlation  matrix  to  converge  quickly  to  the 
DOA  vector  in  the  M-dimensional  reduced  interference  beamspace  data. 

Now,  a  preprocessor  based  on  the  power  minimization  principle  can  only  cancel  jammers/interferers 
above  the  noise  floor.  Smart  jammers  can  lie  at  or  below  the  noise  floor  and  still  disrupt  subsequent 
tracking  and  timing  calculations  related  to  the  GPS  signals.  Thus,  rather  than  do  simple  co-phasal 
beamforming  with  the  DOA  vector  estimate  added  interference  protection  may  be  achieved  by 
doing  adaptive  beamforming  at  the  outputs  of  the  M  ST-PM  preprocessors. 

To  this  end,  let  Rxa;  be  an  M  x  M  spatial  correlation  matrix  (not  cyclic)  formed  at  the  outputs 
of  the  M  ST-PM  preprocessors  (in  the  M-dimensional  reduced  interference  beamspace).  What  we 
propose  is  to  use  the  MSNWF  to  solve  the  system  of  equations.  Weighting  and 

summing  the  M  ST-PM  outputs  with  the  elements  of  thus  computed,  as  shown  in  Figure  9, 
represents  a  per  satellite  combining  that  works  to  suppress  jammers/interferers  at  or  just  below  the 
noise  floof  while  maintaining  an  SWNR  gain  of  M  as  best  as  possible. 

3.7  Power  Minimization  Space-Time  Preprocessor  Followed  by  Per  Satellite  MMSE 
Receiver  Based  on  Known  C/A  Code 

As  discussed  in  Section  3.1,  the  fourth  AJ  processor  presented  incorporates  the  second  and  third 
presented  AJ  processors  as  the  first  two  layers  of  processing,  and  adds  a  third  layer  of  processing. 
The  first  layer  of  processing  entails  multiple  Space-Time  Power  Minimization  (ST-PM)  preprocessors 
as  presented  in  Section  refSTPMBank,  working  to  cancel  jammers/interferences  above  the  noise  floor. 
The  second  layer  of  processing  is  the  per  satellite  combining  of  these  multiple  Space-Time  Power 
Minimization  (ST-PM)  preprocessor  outputs  based  on  an  estimate  of  the  DOA  vector.  This  works 
to  cancel  jammers/interferers  near  the  noise  floor.  The  third  layer  of  processing  entails  a  Minimum 
Mean  Square  Error  (MMSE)  receiver  applied  to  the  output  of  module  2  based  on  the  known  C/A 
Gold  Code.  This  works  to  cancel  jammers/interferers  with  respective  power  levels  below  the  noise 
floor,  including  delay  spoofers  and  multipath  propagation. 

The  presented  MMSE  receiver  works  in  code-space  to  effect  cancellation  of  multi-user  access  in¬ 
terference  (MAI)  in  direct-sequence  spread  spectrum  (DS-SS)  systems  where  the  code  for  each  user  is 
repeated  for  each  and  every  transmitted  information  symbol,  as  is  the  case  for  the  C/A  code  signal. 
The  key  concepts  underlying  the  algorithm  are  explained  in  terms  of  how  it  would  be  applied  in  the 
case  of  the  C/A  code  signal.  Let  denote  an  Nq  x  i  vector  containing  the  Nq  =  1023  length  Gold 
Code  for  the  k-th  GPS  satellite.  Let  B^  denote  an  Nq  x  (Nq  -  1)  blocking  matrix  containing  as 
columns  Nq  —  1  vectors,  each  of  which  is  orthogonal  to  .  Each  column  of  B^  is  viewed  as  an  FIR 
(digital)  filter  as  is  c0°). 

Once  synchronization  with  the  GPS  from  the  k- th  satellite  is  achieved,  each  of  the  respective 
outputs  of  the  bank  of  Nq  —  1  parallel  filters  formed  from  the  columns  of  B^  only  contains  interference 
or  jammers,  assuming  one  slides  over  a  full  Gold  Code  duration  at  a  time.  A  weighted  sum  of  these  is 
thus  used  to  estimate  the  interference  contaminating  the  output  of  the  code-matched  filter  c^kk  The 
receiver  adapts  the  weighting  of  the  respective  outputs  of  the  Nq  —  1  parallel  filters  formed  from  the 
columns  of  Bq^  so  as  to  minimize  the  mean  square  value  of  the  difference  between  the  resulting  sum 
and  the  output  of  c 
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3.7.1  Combining  Partial  Despreading  and  the  Multi-Stage  Nested  Wiener  Filter:  A 
Powerful  Innovation 

The  weight  vector  to  be  adapted  is  essentially  as  long  as  the  code  length.  In  order  to  speed  up 
convergence,  and  to  reduce  computational  complexity,  We  have  applied  the  MSNWF  to  this  problem 
by  defining  the  “desired”  signal  as  the  output  of  the  the  code-matched  filter  c^).  This  yielded  a 
reduced-rank  adaptive  algorithm  for  finding  that  weight  vector  which  maximizes  the  post-correlation 
SINR  for  a  given  “user”  knowing  only  the  code  of  that  “user”,  the  C/A  code  in  our  applicatipn  here. 
This  prior  work  primarily  deals  with  commercial  CDMA  systems,  where  the  channel  code  length  is  on 
the  order  of  a  100.  Here  the  Gold  Code  lenth  is  1023,  an  order  of  magnitude  higher.  Ostensibly,  this 
translates  into  a  very  large  computational  burden. 

However,  as  pointed  out  in  Section  3.4.1,  there  are  relatively  large  Doppler  shifts  due  to  satellite 
motion.  Thus,  to  reduce  the  computational  complexity  AND  counter  the  deleterious  effects  of  Doppler, 
we  propose  to  once  again  employ  the  method  of  partial  despreading.  The  length  1023  Gold  Code  for 
a  given  GPS  satellite  is  decomposed  into  L  nonoverlapping  sub-blocks  of  length  Q ,  with  Q  chosen  as 
large  as  possible  under  the  constraint  that  the  change  in  phase  over  the  sub-block  due  to  Doppler  be 
relatively  negligible.  Since  Q  can  still  be  quite  large,  L  =  16  blocks  of  length  Q  =  64  is  a  reasonable 
starting  point,  we  propose  to  apply  the  MSNWF  to  the  output  of  each  sub-block.  This  combination 
of  partial  despreading  and  the  MSNWF  is  a  truly  powerful  innovation  that  has  not  heretofore  been 
developed  for  any  CDMA  application. 

3.8  Polarization  Diverse  Arrays  and  Compact  Arrays 

As  noted  previously,  the  formation  of  a  spatial  null  towards  an  interferer  may  drop  the  gain  in  the 
direction  of  a  “nearby”  (in  angle)  GPS  satellite  so  low  that  it  is  rendered  useless.  Due  to  the  small 
array  aperture  in  the  GPS  application,  an  interferer  and  a  GPS  satellite  may  be  widely-spaced  in 
terms  of  physical  angle  but  closely-spaced  in  terms  of  beamwidths.  Polarization  diversity  amongst  the 
antenna  elements  is  a  way  to  combat  this  problem. 

For  example,  consider  the  compact  array  pictured  in  Figure  8  (b)  used  for  direction  finding  (DF ) 
and  beamforming  at  HF  for  electronic  warfare  applications.  This  compact  array  consists  of  three, 
orthogonal  loops  and  three  orthogonal  dipoles,  all  co-located  but  having  mutually  orthogonal  orienta¬ 
tions.  Dr.  Zoltowski  has  developed  algorithms  for  DF  and  beamforming  which  exploit  the  polarization 
diversity  inherent  in  this  compact  array  [Zo6,  Zo7].  For  example,  it  has  been  demonstrated  that  two 
signals  coming  from  exactly  the  same  radial  direction  but  with  different  polarization  states  may  be 
extracted  individually! 

Again,  though,  the  compact  array  pictured  in  Figure  8  (b)  was  developed  for  operation  at  HF.  It 
was  presented  here  solely  for  the  purpose  of  emphasizing  the  importance  of  polarization  diversity.  For 
the  GPS  application,  effecting  polarization  diversity  amongst  the  antenna  elements  may  be  as  simple 
as  having  adjacent  antenna  elements  oppositely  polarized,  clock-wise  circular  followed  by  counter 
clock-wise  circular  in  sequence.  Compact  nonlinear  antenna  arrays  and  alternative  realizations  of 
polarization  diversity  have  been  investigated  as  part  of  this  project. 

In  the  GPS  application,  the  beamforming  weights  that  yield  maximum  SINR  for  a  given  GPS 
satellite  have  a  complex  dependence  on  on  the  polarization  states  and  azimuth/elevation  angles  of  both 
that  satellite’s  signal  and  the  interfering  signals.  This  dependence  is  further  complicated  by  mutual 
coupling  effects.  However,  all  of  the  algorithms  presented  for  this  project  adapt  either  based  on  the 
power  minimization  principle  or  the  the  known  code  of  a  given  GPS  satellite.  The  optimal  weights 
are  achieved  automatically  without  the  need  for  estimating  polarization  states  or  azimuth/elevation 
angles. 
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Signal  and  Array  Processing,  SSAP  2000,  Pocono  Manor,  PA,  14-16  August  2000,  pp.  514-518. 

25.  Wilbur  Myrick  and  Michael  D.  Zoltowski,  “Low-Sample  Performance  of  Reduced-Rank  Power 
Minimization  Based  Jammer  Suppression  for  GPS,”  IEEE  Sixth  International  Symposium  on 
Spread  Spectrum  Techniques  &  Applications  (ISSSTA  2000),  Parsippany,  NJ,  6-8  September 
2000,  pp.  93-97. 

26.  W.L.  Myrick,  M.D.  Zoltowski,  and  J.S.  Goldstein,  ’’Adaptive  Anti-Jam  Reduced-Rank  Space- 
Time  Preprocessor  Algorithms  for  GPS,”  Institute  of  Navigation  (ION)  Conference,  Salt  Lake 
City,  Utah,  17-20  Sept.  2000. 

6  INTERACTIONS/TRANSITIONS: 

6.1  A.  Participation/presentations  at  meetings,  conferences,  seminars,  etc 

6.1.1  Invited  Seminars. 

•  “Space-Time  Signal  Processing  for  Wireless  Communications:  Equalization  and  Interference 
Cancellation,”  ECE  Dept.  Colloquium  at  University  of  Minnesota  (televised),  (invited  by  Prof. 
Ahmed  Tewfik)  12  February  1998. 

•  “Smart  Antennas  for  Wireless  Communications,”  (invited  by  Prof.  Mani  Venkata)  Jerry  J unkins 
Endowed  Chair  Seminar  Series,  Department  of  Electrical  and  Computer  Engineering,  Iowa  State 
University,  Ames,  IA,  13  February  1996. 

•  “Semi-Blind  Space-Time  Equalization  for  High-Speed  Linear/Nonlinear  Modulations  Based  on 
Subbanding  and  the  SIMO  Cross-Relation”  (invited  by  Dr.  Timothy  Thomas)  Communications 
Systems  Research  Lab,  Motorola,  Schaumburg,  IL,  18  December  1998. 

•  “Overview  of  Physical  Layer  Standards  and.  Related  Space-Time  Signal  Processing  Research  at 
Purdue  University  for  Third  Generation  Wireless  Communications,”  ECE  Dept.  Colloquium  at 
University  of  Virginia  (Eminent  Speaker  Series),  (invited  by  Prof.  Nikos  Siridopolous)  12  May 
1999. 

•  “Overview  of  Physical  Layer  Standards  and  Related  Space-Time  Signal  Processing  Research  at 
Purdue  University  for  Third  Generation  Wireless  Communications,”  ECE  Dept.  Colloquium  at 
University  of  Illinois  at  Chicago,  (invited  by  Prof.  Arye  Nehorai),  2  June  1999. 

•  “Dual-Channel  Zero  Forcing  Equalization  for  3G  CDMA  Forward  Link  to  Restore  Orthogonality 
of  OVSF  Channel  Codes,”  NCE  Cullimore  Memoral  Distinguished  Lecture  Series  at  New  Jersey 
Institute  of  Technology,  Newark,  NJ,  (invited  by  Prof.  Yeshkel  Bar-Ness),  1  November  1999. 

•  “Space-Time  RAKE  Receiver/Zero  Forcing  Equalizer  Hybrids  for  3G  CDMA  Forward  Link  Ex¬ 
ploiting  Orthogonality  of  OVSF  Channel  Codes,”  ECE  Dept.  Colloquium  at  Nothwestern  Uni¬ 
versity,  Evanston,  IL,  (invited  by  Prof.  Michael  Honig),  12  November  1999. 

6.2  B.  Consultative  and  advisory  functions  to  other  laboratories  and  agencies, 
especially  Air  Force  and  other  DoD  laboratories 

•  Served  on  review  panel  for  National  Science  Foundation’s  Wireless  Information  Technology  and 
Networks  Proposal  Review  Panel,  . CISE  Directorate,  CCR,  Signal  Processing  Systems  Program, 
Arlington,  VA,  21-22  June  1999. 
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6.3  C.  Transitions 

None  to  report  at  this  time.  We  will  continue  our  dialogue  with  Dr.  James  Tsui  at  Wright  Laboratories 
on  the  Wright-Patterson  Air  Force  Base  in  Dayton,  OH  to  work  towards  a  transition  of  this  research. 
Dr.  Tsui  is  currently  overseeing  the  development  of  an  experimental  multichannel  receiver  for  GPS 
at  Wright  Labs.  The  initial  goal  of  this  interaction  is  to  test  the  GPS  anti-jam  algorithms  developed 
through  this  effort  with  experimental  data  supplied  by  Dr.  Tsui.  We  will  also  continue  our  dialogue 
with  Eddie  Gibbs  at  Eglin  Air  Force  Base  to  effect  transitions  of  this  research  as  well.  We  are  also 
working  to  effect  transitions  to  the  GPS  JPO  Office  through  SAIC’s  GPS  group  in  Los  Angeles. 

7  NEW  DISCOVERIES,  INVENTIONS,  OR  PATENT  DISCLO¬ 
SURES 

No  patent  disclosures  to  report.  New  discoveries  and  inventions  summarized  previously. 

8  HONORS/AWARDS 

•  Recipient  along  with  Der-Feng  Tseng  of  “The  Fred  Ellersick  MILCOM  Award  for  Best  Paper  in 
the  Unclassified  Technical  Program”  at  the  IEEE  Military  Communications  (MILCOM  ’98)  Con¬ 
ference  held  in  Bedford,  MA  19-21  October  1998.  The  paper  was  entitled  “Blind  Multichannel 
Identification  for  High-Speed  TDMA”  and  the  plaque  was  awarded  by  the  IEEE  Communica¬ 
tions  Society  and  the  MILCOM  Conference  Board  at  the  Chairman’s  Banquet  held  during  the 
conference.  I  was  the  lead  and  primary  author  of  the  paper  and  the  presenter  as  well.  MILCOM 
’98  is  a  major  IEEE  sponsored  conference  in  Communications  and  there  was  only  one  Paper 
Award  out  of  roughly  230  unclassified  papers. 

•  I  was  elected  Fellow  of  IEEE,  effective  1  January  1999  for  “Contributions  to  the  theory  of  antenna 
array  signal  processing  and  two-dimensional  direction-of-arrival  estimation”.  Photo  and  Fellow 
citation  to  be  included  in  April  1999  issue  of  IEEE  Communications  Magazine. 

•  I  was  elected  Secretary  of  the  IEEE  Signal  Processing  Society,  elected  16  May  1998  by  the  Board 
of  Governors,  for  a  three  year  term  starting  1  January  1999.  As  Secretary,  I  automatically 
became  a  member  of  the  Executive  Committee  of  IEEE  Signal  Processing  Society. 

•  I  was  elected  to  membership  in  Signal  Processing  for  Communications  Technical  Committee 
(SSAP  TC)  of  the  IEEE  Signal  Processing  Society,  elected  for  a  three  year  term  starting  1 
January  1999. 

•  I  continued  to  serve  as  an  elected  Member-at-Large  of  the  Board  of  Governors  of  the  IEEE  Signal 
Processing  Society.  My  three  year  term  from  1  January  1998  and  ends  31  December  2000. 
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